Signatures of Emerging Subsurface Structures in Acoustic Power Maps

We show that under certain conditions, subsurface structures in the solar interior can alter the average acoustic power observed at the photosphere above them. By using numerical simulations of wave propagation, we show that this effect is large enough for it to be potentially used for detecting emerging active regions before they appear on the surface. In our simulations, simplified subsurface structures are modeled as regions with enhanced or reduced acoustic wave speed. We investigate the dependence of the acoustic power above a subsurface region on the sign, depth, and strength of the wave speed perturbation. Observations from the Solar and Heliospheric Observatory/Michelson Doppler Imager (SOHO/MDI) prior and during the emergence of NOAA active region 10488 are used to test the use of acoustic power as a potential precursor of magnetic flux emergence.Comment: 7 pages, 5 figures, accepted for publication in Solar Physics on 21 March 201

Helioseismic Holography of an Artificial Submerged Sound Speed Perturbation and Implications for the Detection of Pre-Emergence Signatures of Active Regions

We use a publicly available numerical wave-propagation simulation of Hartlep et al. 2011 to test the ability of helioseismic holography to detect signatures of a compact, fully submerged, 5% sound-speed perturbation placed at a depth of 50 Mm within a solar model. We find that helioseismic holography as employed in a nominal "lateral-vantage" or "deep-focus" geometry employing quadrants of an annular pupil is capable of detecting and characterizing the perturbation. A number of tests of the methodology, including the use of a plane-parallel approximation, the definition of travel-time shifts, the use of different phase-speed filters, and changes to the pupils, are also performed. It is found that travel-time shifts made using Gabor-wavelet fitting are essentially identical to those derived from the phase of the Fourier transform of the cross-covariance functions. The errors in travel-time shifts caused by the plane-parallel approximation can be minimized to less than a second for the depths and fields of view considered here. Based on the measured strength of the mean travel-time signal of the perturbation, no substantial improvement in sensitivity is produced by varying the analysis procedure from the nominal methodology in conformance with expectations. The measured travel-time shifts are essentially unchanged by varying the profile of the phase-speed filter or omitting the filter entirely. The method remains maximally sensitive when applied with pupils that are wide quadrants, as opposed to narrower quadrants or with pupils composed of smaller arcs. We discuss the significance of these results for the recent controversy regarding suspected pre-emergence signatures of active regions

Numerical MHD Simulations of Solar Magnetoconvection and Oscillations in Inclined Magnetic Field Regions

The sunspot penumbra is a transition zone between the strong vertical magnetic field area (sunspot umbra) and the quiet Sun. The penumbra has a fine filamentary structure that is characterized by magnetic field lines inclined toward the surface. Numerical simulations of solar convection in inclined magnetic field regions have provided an explanation of the filamentary structure and the Evershed outflow in the penumbra. In this paper, we use radiative MHD simulations to investigate the influence of the magnetic field inclination on the power spectrum of vertical velocity oscillations. The results reveal a strong shift of the resonance mode peaks to higher frequencies in the case of a highly inclined magnetic field. The frequency shift for the inclined field is significantly greater than that in vertical field regions of similar strength. This is consistent with the behavior of fast MHD waves.Comment: 9 pages, 6 figures, Solar Physics (in press

Constructing and Characterising Solar Structure Models for Computational Helioseismology

In this paper, we construct background solar models that are stable against convection, by modifying the vertical pressure gradient of Model S (Christensen-Dalsgaard et al., 1996, Science, 272, 1286) relinquishing hydrostatic equilibrium. However, the stabilisation affects the eigenmodes that we wish to remain as close to Model S as possible. In a bid to recover the Model S eigenmodes, we choose to make additional corrections to the sound speed of Model S before stabilisation. No stabilised model can be perfectly solar-like, so we present three stabilised models with slightly different eigenmodes. The models are appropriate to study the f and p1 to p4 modes with spherical harmonic degrees in the range from 400 to 900. Background model CSM has a modified pressure gradient for stabilisation and has eigenfrequencies within 2% of Model S. Model CSM_A has an additional 10% increase in sound speed in the top 1 Mm resulting in eigenfrequencies within 2% of Model S and eigenfunctions that are, in comparison with CSM, closest to those of Model S. Model CSM_B has a 3% decrease in sound speed in the top 5 Mm resulting in eigenfrequencies within 1% of Model S and eigenfunctions that are only marginally adversely affected. These models are useful to study the interaction of solar waves with embedded three-dimensional heterogeneities, such as convective flows and model sunspots. We have also calculated the response of the stabilised models to excitation by random near-surface sources, using simulations of the propagation of linear waves. We find that the simulated power spectra of wave motion are in good agreement with an observed SOHO/MDI power spectrum. Overall, our convectively stabilised background models provide a good basis for quantitative numerical local helioseismology. The models are available for download from http://www.mps.mpg.de/projects/seismo/NA4/.Comment: 35 pages, 23 figures Changed title Updated Figure 1

Newly identified properties of surface acoustic power

The cause of enhanced acoustic power surrounding active regions, the acoustic halo, is not as yet understood. We explore the properties of the enhanced acoustic power observed near disk center from 21 to 27 January 2002, including AR 9787. We find that (i) there exists a strong correlation of the enhanced high frequency power with magnetic-field inclination, with greater power in more horizontal fields, (ii) the frequency of the maximum enhancement increases along with magnetic field strength, and (iii) the oscillations contributing to the halos show modal ridges which are shifted to higher wavenumber at constant frequency in comparison to the ridges of modes in the quiet-Sun.Comment: 16 pages, 10 figures, submitted to solar physic

Determining Absorption, Emissivity Reduction, and Local Suppression Coefficients inside Sunspots

The power of solar acoustic waves is reduced inside sunspots mainly due to absorption, emissivity reduction, and local suppression. The coefficients of these power-reduction mechanisms can be determined by comparing time-distance cross-covariances obtained from sunspots and from the quiet Sun. By analyzing 47 active regions observed by SOHO/MDI without using signal filters, we have determined the coefficients of surface absorption, deep absorption, emissivity reduction, and local suppression. The dissipation in the quiet Sun is derived as well. All of the cross-covariances are width corrected to offset the effect of dispersion. We find that absorption is the dominant mechanism of the power deficit in sunspots for short travel distances, but gradually drops to zero at travel distances longer than about 6 degrees. The absorption in sunspot interiors is also significant. The emissivity-reduction coefficient ranges from about 0.44 to 1.00 within the umbra and 0.29 to 0.72 in the sunspot, and accounts for only about 21.5% of the umbra's and 16.5% of the sunspot's total power reduction. Local suppression is nearly constant as a function of travel distance with values of 0.80 and 0.665 for umbrae and whole sunspots respectively, and is the major cause of the power deficit at large travel distances.Comment: 14 pages, 21 Figure

Time--Distance Helioseismology Data Analysis Pipeline for Helioseismic and Magnetic Imager onboard Solar Dynamics Observatory (SDO/HMI) and Its Initial Results

The Helioseismic and Magnetic Imager onboard the Solar Dynamics Observatory (SDO/HMI) provides continuous full-disk observations of solar oscillations. We develop a data-analysis pipeline based on the time-distance helioseismology method to measure acoustic travel times using HMI Doppler-shift observations, and infer solar interior properties by inverting these measurements. The pipeline is used for routine production of near-real-time full-disk maps of subsurface wave-speed perturbations and horizontal flow velocities for depths ranging from 0 to 20 Mm, every eight hours. In addition, Carrington synoptic maps for the subsurface properties are made from these full-disk maps. The pipeline can also be used for selected target areas and time periods. We explain details of the pipeline organization and procedures, including processing of the HMI Doppler observations, measurements of the travel times, inversions, and constructions of the full-disk and synoptic maps. Some initial results from the pipeline, including full-disk flow maps, sunspot subsurface flow fields, and the interior rotation and meridional flow speeds, are presented.Comment: Accepted by Solar Physics topical issue 'Solar Dynamics Observatory

Modeling the Subsurface Structure of Sunspots

While sunspots are easily observed at the solar surface, determining their subsurface structure is not trivial. There are two main hypotheses for the subsurface structure of sunspots: the monolithic model and the cluster model. Local helioseismology is the only means by which we can investigate subphotospheric structure. However, as current linear inversion techniques do not yet allow helioseismology to probe the internal structure with sufficient confidence to distinguish between the monolith and cluster models, the development of physically realistic sunspot models are a priority for helioseismologists. This is because they are not only important indicators of the variety of physical effects that may influence helioseismic inferences in active regions, but they also enable detailed assessments of the validity of helioseismic interpretations through numerical forward modeling. In this paper, we provide a critical review of the existing sunspot models and an overview of numerical methods employed to model wave propagation through model sunspots. We then carry out an helioseismic analysis of the sunspot in Active Region 9787 and address the serious inconsistencies uncovered by \citeauthor{gizonetal2009}~(\citeyear{gizonetal2009,gizonetal2009a}). We find that this sunspot is most probably associated with a shallow, positive wave-speed perturbation (unlike the traditional two-layer model) and that travel-time measurements are consistent with a horizontal outflow in the surrounding moat.Comment: 73 pages, 19 figures, accepted by Solar Physic

Lifetimes of High-Degree p Modes in the Quiet and Active Sun

We study variations of the lifetimes of high-degree solar p-modes in the quiet and active Sun with the solar activity cycle. The lifetimes in the degree range 300 - 600 and frequency 2.5 - 4.5 mHz were computed from SOHO/MDI data in an area including active regions and quiet Sun using the time-distance technique. We applied our analysis to the data in four different phases of solar activity: in 1996 (at minimum), 1998 (rising phase), 2000 (at maximum) and 2003 (declining phase). The results from the area with active regions show that the lifetime decreases as activity increases. The maximal lifetime variations are between solar minimum in 1996 and maximum in 2000; the relative variation averaged over all mode degree values and frequencies is a decrease of about 13%. The lifetime reductions relative to 1996 are about 7% in 1998 and about 10% in 2003. The lifetime computed in the quiet region still decreases with solar activity although the decrease is smaller. On average, relative to 1996, the lifetime decrease is about 4% in 1998, 10% in 2000 and 8% in 2003. Thus, measured lifetime increases when regions of high magnetic activity are avoided. Moreover, the lifetime computed in quiet regions also shows variations with activity cycle.Comment: 13 pages, 5 figures; Accepted for publication in Solar Physic

Local Helioseismology of Sunspots: Current Status and Perspectives (Invited Review)

Mechanisms of the formation and stability of sunspots are among the longest-standing and intriguing puzzles of solar physics and astrophysics. Sunspots are controlled by subsurface dynamics hidden from direct observations. Recently, substantial progress in our understanding of the physics of the turbulent magnetized plasma in strong-field regions has been made by using numerical simulations and local helioseismology. Both the simulations and helioseismic measurements are extremely challenging, but it becomes clear that the key to understanding the enigma of sunspots is a synergy between models and observations. Recent observations and radiative MHD numerical models have provided a convincing explanation to the Evershed flows in sunspot penumbrae. Also, they lead to the understanding of sunspots as self-organized magnetic structures in the turbulent plasma of the upper convection zone, which are maintained by a large-scale dynamics. Local helioseismic diagnostics of sunspots still have many uncertainties, some of which are discussed in this review. However, there have been significant achievements in resolving these uncertainties, verifying the basic results by new high-resolution observations, testing the helioseismic techniques by numerical simulations, and comparing results obtained by different methods. For instance, a recent analysis of helioseismology data from the Hinode space mission has successfully resolved several uncertainties and concerns (such as the inclined-field and phase-speed filtering effects) that might affect the inferences of the subsurface wave-speed structure of sunspots and the flow pattern. It becomes clear that for the understanding of the phenomenon of sunspots it is important to further improve the helioseismology methods and investigate the whole life cycle of active regions, from magnetic-flux emergence to dissipation.Comment: 34 pages, 18 figures, submitted to Solar Physic